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We develop the idea of selectively manipulating the condensate in a trapped Bose-condensed gas, without perturbing 
the thermal cloud. The idea is based on the possibility to modify the mean field interaction between atoms 
(scattering length) by nearly resonant incident light or by spatially uniform change of the trapping magnetic field. 
For the gas in the Thomas-Fermi regime we find analytical scaling solutions for the condensate wavefunction 
evolving under arbitrary variations of the scattering length a. The change of a from positive to negative induces a 
global collapse of the condensate, and the final stages of the collapse will be governed by intrinsic decay processes. 

PACS numbers: 34.20.Cf, 03.75.Fi 



The discovery of Bose-Einstein condensation (BEC) in 
trapped clouds of ultra-cold alkali atoms opened 
unique possibilities to investigate collective many-body 
effects in dilute gases. Isolation of a trapped gas from 
the environment, provided by a wall-less magnetic con- 
finement, makes this object attractive for studying fun- 
damental problems of the physics of many-body quantum 
systems, such as relaxation and the loss of coherence in 
the evolving macroscopic quantum state (condensate). 
In ongoing experiments the condensate is set into 

motion (for example, undergoes oscillations) by varying 
the confining field. Scaling theory of coherent evolution 
(without damping and relaxation) of a Bose condensate 
in harmonic traps under arbitrary frequency variations 
has been developed in PJTo|. It is important, however, 
that variations of the confining potential also cause the 
evolution of the thermal component of a trapped gas. 
Moreover, in the hydrodynamic regime the thermal cloud 
in many aspects evolves similarly to the condensate ||. 
For example, the asymmetry of free expansion and eigen- 
frequencies of small oscillations are almost the same. 

In this paper we develop the idea of selectively manip- 
ulating the condensate, without perturbing the thermal 
cloud. This is especially important for studying the in- 
teraction between the condensate and the thermal com- 
ponent. The idea is based on the possibility to mod- 
ify the mean field interaction between atoms. At low 
temperatures the latter is proportional to the scattering 
length a which, as found in [pd| , can be changed under 
the influence of red-detuned nearly resonant light. An- 
other option to modify the scattering length relies on the 
magnetic field dependence of a, predicted in |l2), and as- 
sumes spatially uniform variations of the trapping field, 
without changing the trap frequencies. Since the shape 
of the condensate wavefunction is predetermined by the 
interaction between particles, the change of a will cause 
the evolution of the condensate at constant frequencies 
of the confining potential. Remarkably, at temperatures 
T nil (n is the gas density, U — 4:Trh 2 a/m, and m 
the atom mass) only the condensate evolution will be 
pronounced, which resembles the picture of the fourth 



sound in superfluid helium. In this temperature range 
the shape of the thermal cloud is mainly determined by 
the trapping potential and temperature. The perturba- 
tion of the thermal component will be small compared to 
that of the condensate, at least as nU /T. 

The change of a in the field of nearly resonant red- 
detuned light is provided by virtual radiative transitions 
of a pair of atoms to a bound electronically excited molec- 
ular state. Since the resonance dipole interaction be- 
tween atoms in the excited state is much stronger than 
the interaction in the ground state, already for moderate 
light intensities the effective interatomic interaction and, 
hence, the scattering length can be significantly changed. 
It is even possible to switch the sign of a |ll]] . 

The magnetic field dependence of the scattering length, 
found in |T^ ], has a resonance structure. Therefore, spa- 
tially uniform variations of the field can also change a 
substantially and reverse its sign. 

Below we find analytical scaling solutions for the con- 
densate wavefunction evolving under arbitrary variations 
of the scattering length. We consider the most inter- 
esting situation, where the mean field interaction be- 
tween atoms greatly exceeds the level spacing in the trap 
(Thomas- Fermi regime). The case of switching a from 
positive to negative and inducing a global collapse of the 
condensate is specially analyzed. 

Let us consider the evolution of a condensate in a har- 
monic potential V(r) = mJ2i UJ i r i /%> with constant fre- 
quencies u)i, induced by variations of the mean field inter- 
action between particles. The equation for the conden- 
sate wavefunction \&o(r, t), with above-condensate parti- 
cles neglected, can be represented in the form 
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dt 



.|_A* + ^£u;M*o + tf(t)|*o| 3 *o. (1) 



Here U(t) — 4:Trh 2 a(t)/m, and a(t) is the time-dependent 
scattering length. Along the lines of the general method 
developed in Q , we introduce d scaling parameters bi (t) , 
where d is the dimension of the system. Turning to 
new coordinate and time variables pi = ri/bi(t), r(t) = 



1 



fdt'f(t')/V(t'), with 7 (t) = U(t)/U and U the initial 
value of U(t), we search for the solution in the form 

* (r, t) = V- 1 / a (t)xo(p, r(f)) exp (i*(r, i)), (2) 

where the dimensionless volume V(i) = ]X bi(t). Substi- 
tuting Eq.(|^) into Eq.(^j) we require the cancellation of 
the VpXo terms. This leads to the relation for the phase: 

$(r,t) = (m/270l>?[M*)/M*)]- (3) 

The scaling parameters bi(t) will be chosen such that they 
are governed by the equations 

S i + w?6 i =w i 2 7(t)/6 i V(t), (4) 

with initial conditions 6^(0) = 1, &i(0) = 0. Then, wc 
arrive at the equation of motion 

^ =- ph>+UolA - (5) 

In the case of initially repulsive (a(0) > 0) and strong 
interaction between particles (Thomas- Fermi regime) the 
initial chemical potential fio = HoUq 3> TiuJi (no is the 
maximum density in the initial static condensate), and 
the kinetic energy term in Eq.(^|) is initially small com- 
pared to the non-linear interaction term . In the course 
of evolution the ratio of the kinetic to interaction term 
scales as e(t) — Ylifa^i/ M(A(*)) 2 V(i)/7(i). Assuming 
that the condition e(t) <gc 1 is satisfied at any t, the ki- 
netic energy term can be omitted. In the case of a fast 
change of the interparticle interaction this, in particular, 
requires the inequality | 7 | 3> (fum / fJ-o) 2 at t close to 0. 
Then in the variables pi, r Eq.(j|) is reduced to that for 
the static case, with initial interaction between particles. 
The solution has the well known form (see fl3] , |l4|| ) : 

Xo(p,r(t)) = ^^-^^(=^) (6) 

in the spatial region where the argument of the square 
root is positive and zero otherwise. Thus, the conden- 
sate evolution under arbitrary variations of the scatter- 
ing length is described by a universal scaling solution for 
^o(r,t), following from Eqs. (||) and (0). Actually, the 
problem is reduced to the solution of classical equations 
(||) for the scaling parameters bi(i) (cf. ||). 

We should emphasize that the results for d = 2, ob- 
tained below, can be applied equally well to the radial 
evolution of in very long samples where the axial fre- 
quency is much smaller than the radial one. This follows 
from the fact that then to a first approximation the de- 
pendence of on the axial coordinate can be omitted. 

In the limiting case of adiabatically slow variations of 
the scattering length, that is on a time scale tq 3> u^ 1 , 
the solution of Eqs. (El) leads to simple relations 



bi(t) = [7(i)] 1/(2+d) = &(*); V(t) = [ 1 (t)] d ^ +d \ (7) 

Assuming that at t^$>To the scattering length acquires a 
constant value a\ (U{t) = U\, ~f(t) = Ui/Uo=Ji), for the 
phase in Eq.(^) we obtain Hor(t) = [lit, where the quan- 
tity /ii = /io7 2 ^ 2+<i ''is equal to the chemical potential of 
the gas with the interparticle interaction U\ . Then, from 
Eqs. (§) and @ one can sec that the initial condensate is 
adiabatically transformed to the static condensate with 
the wavefunction corresponding to the interaction U\. 

For an abrupt change of the interparticle interaction 
from Uq to Ui scaling equations (Bj) take the form 

k + ^bi = Ljf 7l /biV(t). (8) 

In fact, for 7i = Ui/Uo > these equations are similar to 
those for the scaling parameters of the condensate evo- 
lution after an abrupt change of the trapping potential, 
with initial frequencies being uji^/ji and final frequen- 
cies LUi (see Q). The solution of Eqs. (||) gives oscillating 
functions bi (t) which ensure undamped oscillations of the 
condensate wavefunction. For 71 close to 1 the frequency 
spectrum of the condensate oscillations coincides with 
the set of eigenfrequencies of small shape oscillations of 
the initial condensate, found for cylindrically symmetric 
traps in the JIL A [H and MIT Q| experiments and cal- 
culated in Il5|jl^ , |l0[ p| ■ It is important that an abrupt 
change of the scattering length is equivalent to equal rel- 
ative change of all frequencies. This reduces the possibil- 
ity of stochastization of the condensate evolution in an 
anisotropic harmonic potential, compared to the case of 
independent change of the frequencies. 

The structure of the condensate oscillations caused by 
the change of the scattering length can be easily analyzed 
in the case of isotripic trapping configuration. For d = 2 
Eq.(||) has analytical solution 

bj{t) = b 2 {t) = [1 + 7l + (1 - 7l ) cos 2wt]/2. (9) 

Once the scattering length is decreased (71 < 1) there 
will be "compression oscillations": Compression of the 
condensate will be followed by its expansion to the ini- 
tial shape. An increase of a (71 > 1) induces "expansion 
oscillations" , i.e. the condensate is expanding and then 
compressing to the initial shape. In both cases the scal- 
ing parameter b varies from 1 to a/7i- In the 3-d case 
the condensate oscillations are anharmonic, with a char- 
acteristic period somewhat smaller than that for d = 2. 
The parameter b is varying from 1 to -y/2 7 i/3. 

The dynamics of the system, induced by the change of 
the interparticle interaction, is drastically different from 
that caused by variations of the trap frequencies. In the 
former case the condensate oscillates, whereas the ther- 
mal component is only weakly perturbed and practically 
remains static. In the latter case the thermal cloud is 
fully involved in the oscillatory evolution process. Pref- 
erential evolution of the condensate under the change of 
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the scattering length resembles the fourth sound in super- 
fluid helium, where the superfluid part of the density os- 
cillates on the background of a static normal component. 
As well as in helium, the analysis of dynamic properties 
of a trapped gas in these conditions can be an effective 
method of identifying and studying BEC. 

The change of the interaction between particles can 
lead to an interesting phenomenon, a global collapse of 
the condensate as a whole. Let us assume that initially 
the scattering length a > 0, and the system is in a sta- 
tionary state. The sign of a can be switched to negative 
in the field of nearly resonant red-detuned light |ll]] or 
by spatially uniform variations of the confining magnetic 
field (see jl2)). Then, for an abrupt change of a the solu- 
tion of scaling equations (||) immediately leads to a self- 
compression (collapse) of the condensate. In the d = 2 
isotropic case this directly follows from the exact solu- 
tion (||): For 71 < the scaling parameter b decreases 
and reaches zero at t = t*. where 



t* = u) 1 arcsin 
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(10) 



In the 3-d isotropic case t* is close to that from Eq.([To 
In the vicinity of t*, where wAt-Cmin^ip/f^il -1 ^] (At: 
f* — i), the solution of Eq.(EJ) takes the form 



b 2 (t) 



^,At)(2 + d)/V2d] 4/(2+d) . (11) 



In both cases the compression rate increases with de- 
creasing At: The quantity b{t)/b{t) = 2/[(d + 2)At], and 
the characteristic velocity of the condensate boundary, 
v cx [At]~ d /( d+2 ). In the course of collapse the small pa- 
rameter e(t) remains constant (d — 2) or even decreases 
(d = 3), which justifies the neglect of the kinetic energy 
term in Eq.(|^). The kinetic energy of the system is de- 
termined by the behavior of the phase $(r) (||) which 
drastically increases with decreasing At. A strong rise 
of the kinetic energy in the course of collapse is com- 
pensated by decreasing potential energy (a < 0), which 
ensures the conservation of the total energy. 

The increase of density in the collapsing condensate 
enhances intrinsic inelastic processes, such as three-body 
recombination and spin-dipole relaxation in binary colli- 
sions. This leads to particle losses, since fast atoms and 
molecules produced in the inelastic processes escape from 
the trap. We will confine ourselves to the analysis of the 
influence of inelastic processes on the dynamics of col- 
lapse in the isotropic d = 2 and d = 3 cases. The atom 
loss rate is determined by the relation 



N 



, J d d r\V Q (r, t)| 4 - a r J d d r\^ Q (r, t)| 6 , (12) 



where a s and a r are the rate constants of spin relax- 
ation and three-body recombination, respectively. With 
increasing density in the course of collapse, the ratio of 
the recombination to relaxation rate increases. In ultra- 
cold alkali atom gases a s is in the range 10" 14 — 10~ 16 



cm 3 /s and a r ~ 10~ 28 cm 6 /s. Accordingly, considering 
initial condensate densities Uq ~ 10 13 — 10 14 cm" 3 (cur- 
rently achieved in Rb and Na BEC experiments 0-0), 
we see that already for moderate compression of the con- 
densate three-body recombination dominates over spin 
relaxation. Therefore, the process of spin relaxation is 
omitted in the further analysis. For the same reason, in 
the case of light-induced change of the scattering length 
we omit the process of photoassociation in pair collisions. 

The loss of particles from the condensate influences the 
evolution of the condensate wavefunction. The character 
of the evolution depends on the density of the collaps- 
ing condensate, n(t), and is predetermined by the re- 
lation between the characteristic time of three-body re- 
combination, r r ~ (a r n , ) , and the correlation time 
r c ~ %/n\U\\. The latter is responsible for establishing 
the shape of the condensate wavefunction correspond- 
ing to the instantaneous value of the number of particles 
N(t). There will be two different stages of the evolution. 
The role of three-body recombination becomes especially 
important at densities 

n > = \Ui\/h,a r ^> no, (13) 

where the compression is very strong and r r < r c . Then 
already the shape of |\I , o( r )i)| is determined by the re- 
combination losses (see below). With realistic numbers 
|ai| ~ 10 A, and a r ~ 10 -28 cm 6 /s, we see that the den- 
sity n* will be in the range 10 17 — 10 18 cm -3 and still sat- 
isfies the criterion of weakly interacting gas, n»|a| 3 <C 1. 

In the density range n <C n* the recombination time 
TV 3> t c and the loss of particles occurs in a quasistation- 
ary regime. For ^o(r, t) one can use the scaling solution 
following from Eqs. (0), (||), with instantaneous values of 
N and the chemical potential fj,. Since /i ~ J\f 2 ^ 2+d \ for 
the time dependence of the maximum condensate density 
in the quasistationary stage we obtain 



n(t) = n Q [b(t)}- d [N(t)/N } 2 ^ 2+d \ 



(14) 



where Nq is the initial number of condensate particles. 

Below we assume that |7i| < 1. Then the characteristic 
time of collapse, t» ~ cj" 1 . Since for realistic parameters 
of the system the initial recombination time r r o 3> cj" 1 , 
a major part of the compression occurs without particle 
losses. This is the case even at t rather close to t*, the 
condensate density being determined by Eq. dl4| ) with the 
scaling parameter b from Eq.(pl|) and N « Nq: 

n(t) w n (VM^t)- 2d ^ 2+d l (15) 

Eq. ( |l5|) is no longer valid when the number of lost par- 
ticles becomes comparable with N . The characteristic 
time Ati and the density iij, at which this happens can 
be found from Eq.([l2|), with ^ determined by Eqs. (||), 
(||) and the scaling parameter b(t) given by Eq.(^Tj): 

d d p 



ui 



N =- TT, I5w77T (MO- mc.y/2) 3 w- - 



b 2d (t) 



HwAt) 4d /( 2 + d ) ' 
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Comparing the number of lost atoms with Nq we obtain 

y/\^\uAt L » (VTriT^o)- (2+d)/(3d - 2) ; (16) 

n L ^n (^\LUT r0 ) 2d/{3d - 2) . (17) 



Assuming \/Yfi\u)T r o > 1 we see that rij, ^> no. 

Eq.(p"3|) can be rewritten as n*«no|7i|(wT r o)(Mo/^)) 
and the condition |7i| 3> (Tiw/rto) 2 (see above) ensures the 
inequality riL^n* which justifies the use of the quasis- 
tationary approach in deriving Eqs. (|l5|), ( |l6| ) and (p"7|). 

The condition nj, <C H* means that the collapse contin- 
ues to occur in the quasistationary regime also at times 
At < At L . In this time interval, using Eqs. (pi]), (14), 
©, and Eq.dH) written in the form N « -a r n 2 (t)N(t), 
for the number of particles wc find 



JV(i) a 7V (Ai/A^)( 3d - 2 )/ 4 . 
The condensate density increases as 



i(f) « n y/T r o/At, 



(18) 



(19) 



i.e. slower than in the initial stage of the compression. 

At times very close to t*, i.e., at At « r r * = (a,-nj) _1 , 
the density approaches n*, and the picture of collapse 
changes. If n becomes larger than n*, then the decrease 
of n due to particle losses already dominates over the 
dynamic compression. On the other hand, for n < n* the 
dynamic compression dominates and the density rises. 
Therefore, n* proves to be a critical density value which 
is approximately conserved in the final stage of collapse. 
The decrease of the number of particles is determined by 
the relation N = — iV/r r *, and we obtain 



N(t) « TV* exp {-(t - Q/r r *}, 



(20) 



where N* <C Nq and is determined by Eqs. (|T^) and (16), 
with At « T r *. The characteristic size of the condensate 
decreases as N(t)/n*. 

The above described picture of collapse is somewhat 
idealized. For example, the mean interparticle separation 
at limiting density, n* x ^ 3 , can become comparable with 
the characteristic radius of interatomic interaction. Nev- 
ertheless, we believe that the main results should remain 
unchanged. The scenario of collapse can be observed by 
measuring the recombination losses of particles. 

Considering the global collapse in an anisotropic har- 
monic potential we should return to scaling equations 
(B), with 7i < 0. The rate of the self-compression is the 
highest in the direction corresponding to the largest fre- 
quency. In the case of cylindrical symmetry with the ratio 
of the axial to radial frequency, (3 = lo z /lo p < 1, there will 
be a characteristic time f* at which b p (t*) = 0, but b z (t*) 
remaines finite. In the time interval, where b p (t) <C 1, 
Eq.(||) for b p (t) is reduced to that for the isotropic 2-d 
collapse, with 71 replaced by 7i/6 z (t*). Thus, the global 
collapse will occur along the lines of the above described 



2-d scenario, in combination with a "slow" compression 
in the axial direction. For (3 > 1 the condensate pre- 
dominantly collapses in the axial direction and, due to a 
quasi-one-dimensional character of the self-compression, 
the kinetic energy term in Eq.(^|) increases. Eventually 
the parameter e(t) becomes of order unity, and our scal- 
ing approach breaks down. Such a rise of the kinetic 
energy prevents the further compression and can lead to 
a non-trivial oscillatory character of the global collapse. 

Density fluctuations on the background of the global 
collapse can lead to the appearence of local collapses, 
with linear dimensions of order the healing length. Al- 
though their influence on the global collapse can be re- 
duced by decreasing the ratio uqU /Thj, the study of the 
role of the local collapses requires a separate analysis. 

We acknowledge discussions with J.T.M. Walraven and 
M.W. Reynolds. This work was supported by the Dutch 
Foundation FOM, by NWO (project 047-003.036), by IN- 
TAS, and by the Russian Foundation for Basic Studies. 



[1] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wie- 

man, and E.A. Cornell, Science, 269, 198 (1995). 
[2] C.C. Bradley, C.A. Sackett, J.J. Tolett, and R.G. Hulet, 

Phys. Rev. Lett., 75, 1687 (1995). 
[3] K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van 

Druten, D.S. Durfee, D.M. Kurn, and W. Ketterle, Phys. 

Rev. Lett., 75, 3969 (1995). 
[4] D.S. Jin, J.R. Ensher, M.R. Matthews, C.E. Wieman, 

and E.A. Cornell, Phys. Rev. Lett., 77, 420 (1996). 
[5] M.R. Andrews, M.-O. Mewes, N.J. van Drutten, D.M. 

Kurn, D.S. Durfee, C.G. Townsend, and W. Ketterle, 

Science, 273, 84 (1996). 
[6] M.-O. Mewes, M.R. Andrews, N.J. van Drutten, D.M. 

Kurn, D.S. Durfee, and W. Ketterle, Phys. Rev. Lett., 

77, 416 (1996). 
[7] M.-O. Mewes, M.R. Anderson, N.J. van Drutten, D.M. 

Kurn, D.S. Durfee, C.G. Townsend, and W. Ketterle, 

Phys. Rev. Lett., 77, 988 (1996). 
[8] D.S. Jin, M.R. Matthews, J.R. Ensher, C.E. Wieman, 

and E.A. Cornell, Phys. Rev. Lett., 78 (1997). 
[9] Yu. Kagan, E.L. Surkov, and G.V. Shlyapnikov, Phys. 

Rev. A, 54, R1753 (1996); ibid, 55, R18 (1997). 
[10] Y. Castin and R. Dum, Phys. Rev. Lett., 77, 5315 (1996). 
[11] P.O. Fedichev, Yu. Kagan, G.V. Shlyapnikov, and J.T.M. 

Walraven, Phys. Rev. Lett., 77, 2913 (1996). 
[12] E. Tiesinga, B.J. Verhaar, and H.T.C. Stoof, Phys. Rev. 

A, 47, 4114 (1993). 

[13] V.V. Goldman, I.F. Silvera, and A.J. Leggett, Phys. Rev. 

B, 24, 2870 (1981). 

[14] DA. Huse and ED. Siggia, J. Low Temp. Phys., 46, 137 
(1982). 

[15] S. Stringari, Phys. Rev. Lett., 77, 2360 (1996). 

[16] M. Edwards, PA. Ruprecht, K. Burnett, R.J. Dodd, and 

C, W. Clark, Phys. Rev. Lett., 77, 1671 (1996). 



4 



